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We present the Nf = 2 + 1 clover fermion lattice QCD calculation of the nucleon strangeness form 
factors. We evaluate disconnected insertions using the Z(4) stochastic method, along with unbiased 
subtractions from the hopping parameter expansion. We find that increasing the number of nucleon 
sources for each configuration improves the signal significantly. We obtain G%f{0) = —0.017(25) (07), 
where the first error is statistical, and the second is the uncertainties in and chiral extrapolations. 
This is consistent with experimental values, and has an order of magnitude smaller error. 

PACS numbers: 13.40.-f, 12.38.Gc, 14.20.Dh 



I. INTRODUCTION 

The structure of the nucleon plays an essential role 
in understanding the dynamics of QCD, and experi- 
ments have been providing many intriguing and unex- 
pected results for decades. In particular, the strangeness 
content of the nucleon attracts a great deal of interest 
lately. As the lightest non-valence quark structure, it 
is an ideal probe for the virtual sea quarks in the nu- 
cleon. Extensive experimental/theoretical studies indi- 
cate that the strangeness content varies significantly de- 
pending on the quantum number carried by the ss pair: 
the scalar density is about 10-20% of that of up, down 
quarks, the quark spin is —10 to 0% of the nucleon, and 
the momentum fraction is only a few percent of the nu- 
cleon. Recently, intensive experiments have been car- 
ried out for the electromagnetic form factors by SAM- 
PLE dl, A4 0, HAPPEX 3], and GO d], through parity- 
violating electron scattering (PVES). The global analy- 
ses 0, HQ have produced, e.g., G%{Q'^) = -0.008(16) 
and GljiQ^) = 0.29(21) at = O.lGeV^ [|, but sub- 
stantial errors still exist so that the results are consistent 
with zero. Making tighter constraints on these form fac- 
tors from the theoretical side is one of the challenges in 





FIG. 1: Two different representations for the 3pt function. 
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QCD calculation. Moreover, such constraints, together 
with experimental inputs, can lead to more precise de- 
terminations of various interesting quantities, such as 
the axial form factor G% 0, and the electroweak radia- 
tive corrections including the nucleon anapolc moment, 

Ga UM- 

Unfortunately, the theoretical status of strangeness 
form factors remains quite uncertain. For instance, the 
values for the magnetic moment G\j{0) from different 
analyses vary widely: —0.31(9) in the dispersion relation 
(DR) with pole ansatz ^, —(0.15—0.51) in DR with scat- 
tering kaon clouds [lC|, 0.035 in the quark model ill I. 
and 0.08 — 0.13 in the chiral quark-soliton model [l2|. 
The analyses for the electric form factor are similarly 
ambiguous. Under these circumstances, the most desir- 
able study is the first-principle calculation in QCD, such 
as the lattice simulation. However, the calculation re- 
quires the evaluation of the disconnected insertion (DI) 
(Fig. [1] (right)), which is much more difficult compared 
to the connected insertion (CI) (Fig. [1] (left)) calculation, 
because the straightforward calculation of DI requires all- 
to-all propagators, and is prohibitively expensive. Con- 
sequently, there are only few DI calculations, where the 
all-to-all propagators are stochastically estimated. The 
first calculation was done in the quenched approximation 
with Wilson fermion [H] and gave G|^(0) = -0.36(20), 
and —0.28(10) from the updated calculation jT3|. An- 
other quenched DI calculation with Wilson fermion [Tsj 
obtained G^^(O.lGeV^) = 0.05(6). There are also sev- 
eral indirect estimates using quenched [l^ [l3| or un- 
quenched [11] lattice data for the CI part, and the ex- 
perimental magnetic moments (or electric charge radii) 
for octet baryons as inputs under the assumption of 
isospin symmetry. These estimates obtained G%j{0) = 
-0.046(19) [3 and G%i{0) = -0.066(26) 

In this paper, we provide the first full QCD lattice sim- 
ulation of the direct DI calculation with high statistics. 
The paper is organized as follows. In Sec. |TT1 the for- 
mulation and parameter setup of the lattice calculation 
are presented. Technical details to improve the DI eval- 
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TABLE I: The setup parameters as well as hadron masses. 
Nnoise IS the number of noises in the stochastic estimate and 
Nsrc is the number of nucleon sources for each configuration. 



f^ud 


Nconf 


Nnoise 


Nsrc 


m-n-a rriKa rriNCi 


0.13760 


800 


600 


64 


0.5141(5) 0.5141(5) 1.0859(12) 


0.13800 


810 


600 


82 


0.4302(6) 0.4540(5) 0.9623(16) 


0.13825 


810 


800 


82 


0.3717(7) 0.4141(6) 0.8844(20) 



uation are also given. In Sec. IIIIl we present the lattice 
QCD data, and carefully examine the possible systematic 
uncertainties. Sec. IIVI is devoted to the summary of our 
results. 



II. FORMALISM 

We employ Nf = 2+1 dynamical configurations 
with nonperturbatively 0(0) improved clover fermion 
and RG-improved gauge action generated by CP- 
PACS/JLQCD Collaborations 0. We use /? = 1.83 
and Csw — 1.7610 configurations with the lattice size of 
X T — 16"^ X 32. The lattice spacing was determined 
using X-input or </)-input [19]. We use the averaged scale 
of a^^ — 1.625GeV hereafter, as the uncertainty of the 
scale is negligible compared to the statistical error in our 
results. For the hopping parameters of u,d quarks (K„d) 
and s quark (ks), we use Kud = 0.13825, 0.13800, and 
0.13760, which correspond to = 0.60, 0.70, and 0.84 
GeV, respectively, and Ks = 0.13760 is fixed. We perform 
the calculation only at the dynamical quark mass points, 
with the periodic boundary condition in all directions. 

We calculate the three point function (3pt) IIj^* and 
two point function (2pt) II^p* defined as 

nj^*(p,t2; q,ti; p' ^ p - q,to) ^ ^ 
g+.f.(Si-.-o)(0|T[xjv(x2,t2)JM(2?l,tl)XJv(xo,to)]|0)(l) 



-ip-(x2~xa) , 



tfP'(p,t; to) 



|T[xjv(x,i)Xiv(a'o,io)]|0), (2) 



where xn = ^ahc(j^G^^d\j)uc is the nucleon inter- 
polation field and the insertion is given by the 
point-split conserved current Jp(a; -f- /^/2) = (l/2)x 
{s{x + At)(l + i,,)lJl{x)s{x) - s(x)(l - -i^)U^(x)s(x + . 

Electromagnetic form factors can be obtained using 
■p — 0, p' — —(f kinematics for the forward propagation 
{ti S> ti ^ to) [1^. In this work, we consider the back- 
ward propagation (t2 <C ti <C to) as well, in order to 
increase statistics. The formulas for Sachs electric (mag- 
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2: R\i (upper) and 7?^ (lower) with K„d = 0.13760, 
q'^ = 2-{2ti / LaY , Nsrc ~ 64 (circles) and Nsrc = 4 (triangles, 
with offset for visibility), plotted against the nucleon sink 
time t2. The dashed line is the linear fit where the slope 
corresponds to the form factor. 



netic) form factors G|; {G%,j) are summarized as 



Tr 



Tr[r±-n2pt(±(f,ii; to)] 



Tr 


r± 

^ e 


n2pt(0,ii; to) 


Tr 


J- e 


•n2pt(0,i2; to) 



Rl 



4(rp"oi = rf) 



±G|(Q2), 



pol 



(3) 

(4) 
(5) 



(i±74)/2,r± 



(±i)/2x (1± 



+ q^. The upper (lower) sign 



where {i, j, k} 7^ 4, 
74)757fe and E^; = ^\ 

corresponds to the forward (backward) propagation. 

Furthermore, we consider another kinematics oi p — q, 
p' = 0, where the analogs of Eqs. ([3]), g]), ([5]) hold. We 
find that the results from the latter kinematics have sim- 
ilar size of statistical errors as those from the former, and 
the average of them yields better results. Hereafter, we 
present results from total average of two kinematics and 
forward/backward propagations, unless otherwise noted. 

The calculations of 3pt functions for the strangeness 
current need the evaluation of DI. We use the stochastic 
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TABLE II: Lattice results for the strangeness electromagnetic form factors with the momentum-squared = n ■ (27r/La)^. 





1 


2 3 


4 





GUQ'){xW-^) 
1 2 


3 


4 


0.13760 
0.13800 
0.13825 


-0.96(29) 
-0.76(37) 
-1.09(41) 


-0.61(17) -0.57(20) 
-0.76(24) -0.57(32) 
-1.02(27) -0.67(33) 


-0.20(25) 
0.04(41) 
-0.25(47) 


0.21(21) 
0.15(28) 
-0.04(33) 


0.01(10) 0.14(08) 
0.16(14) 0.14(12) 
-0.16(15) 0.36(13) 


0.22(11) 
0.15(19) 
0.27(20) 


0.10(15) 
0.46(29) 
0.71(31) 



TABLE III; The parameters fitted against behavior. 





GIM 
(xlO~") 


Aa 


x'/dof 


G% 

9e 
(xlO-2) 


iQ') 
xVdof 


0.13760 
0.13800 
0.13825 


-1.7(12) 
-1.4(09) 
-1.9(11) 


0.66(29) 
0.77(40) 
0.80(40) 


0.34(83) 
0.77(126) 
0.38(87) 


1.2(5) 
1.4(6) 
1.9(7) 


0.56(89) 
0.34(70) 
2.68(189) 



method [20|, with Z(4) noises in color, spin and space- 
time indices. We generate independent noises for dif- 
ferent configurations, in order to avoid possible auto- 
correlation. To reduce fluctuations, we use the charge 
conjugation and 75-hermiticity (CH), and parity symme- 
try. For instance, we find that the information for the 
G%i is coded in the product of Re(n^P*) x Re(loop), and 
filtering out the imaginary parts reduces the noises 
[2^ . We also perform unbiased subtractions [1^ to re- 
duce the off-diagonal contaminations to the variance. For 
subtraction operators, we employ those obtained through 
hopping parameter expansion (HPE) for the propagator 
Af-i ,23], ^M-i = j^ + j^{nD)^ + ■■■ where 
D denotes the Wilson-Dirac operator and C the clover 
term. We subtract up to order (k-D)"* term, and observe 
that the statistical error is reduced by a factor of 2. 

In the stochastic method, it is quite expensive to 
achieve a good signal to noise ratio (S/N) just by in- 
creasing Nnoise because S/N improves with yJNnoise- In 
view of this, we use many nucleon point sources iV.srr- in 
the evaluation of the 2pt part for each configuration [24I . 
Since the calculations of the loop part and 2pt part are 
independent of each other, this is expected to be an ef- 
ficient way. In particular, for the Nnoise 2> Ngrc case, 
we observe that S /N improves almost ideally, by a factor 
of yJNsrc- We take Nsrc — 64 for Kud = 0.13760 and 
Nsrc = 82 for Kud = 0.13800,0.13825, where locations 
of sources are taken so that they are separated in 4D- 
volume as much as possible. The calculation parameters 
as well as basic hadron masses are tabulated in Tab. HI 

There are several ways [3, [lH| to extract the matrix 
elements from Eqs. ^ and ^ with various ii, t2 results. 
Among them, it is advocated [H, [H, [2l|, [l^l to take the 
summation over the insertion time ti , symbolically given 
as = T!t\=tl+t, = const, -f x 0% ^, 

E ,M 

where are trivial kinematic factors appearing in 

Eqs. ([4]) and ([5]), and tg is chosen so that the error is 
minimal. We thus obtain G% as the linear slope of 



against t2 ■ In order to achieve ground state satu- 
ration, we use the data only for t2 > 7 [2^ . 

III. RESULTS AND DISCUSSIONS 

We calculate for the five smallest momentum-squared 
points, — n-{2TT / LaY {n ~ 0-4), which correspond to 
\q\ = 0,0.64,0.90,1.1,1.3 GeV. Typical figures for R\.j, 
i?^ are shown in Fig. O One can observe the significant 
S/N improvement by increasing Nsrc- The numerical re- 
sults are given in Tab. |TT1 We note G'|;(0) are consis- 
tent with zero, which serves as a test of the calculations. 
Of particular interest is that, for all n^d simulations, 
G\.j{Q^) is found to be negative with 2-3 a signals for 
low regions. 

In order to determine the magnetic moment, the 
dependence of G\j{Q^) is studied. We employ the dipole 
form in the fit, Gi^(Q2) ^ G|^(0)/(1 + Q'^/K^f, 
where reasonable agreement with lattice data is ob- 
served. For the electric form factor, we employ G'|;(Q^) = 
ah ■ <3V(1 + QVA2)2, considering that G^O) = from 
the vector current conservation. In the practical fit of 
G%{Q'^), however, reliable extraction of the pole mass A 
is impossible because G^^{Q^) data are almost zero within 
error. Therefore, we assume that G%{Q'^) has the same 
pole mass as G\i{Q'^), and perform a one-parameter fit 
for g|;. The obtained parameters are given in Tab. IIIII 
We also test the simultaneous fit of G\,j{Q'^) and G'^^{Q^) 
with three parameters of G|,j(0), A, g'^^, and confirm that 
the results are consistent with the values in Tab. IIIII 

Finally, we perform the chiral extrapolation for the 
fitted parameters. Since our quark masses are rela- 
tively heavy, we consider only the leading dependence 
on rriK, which is obtained by heavy baryon chiral per- 
turbation theory (HBxPT). For the magnetic moment 
G^^(O), we fit Hnearly in terms of rriK [H HJ. For 
the pole mass A, we take that the magnetic mean- 
square radius {r^s)M = ~G^\q2^o = l2Gli{0)/A^ be- 
haves as l/rriK ^24^. For g^^, we use the electric radius 
(^s)b = -6i|f Iq2=o = -65I which has an ln{mK/n) 
behavior [2J|, and we take the scale /i = 1 GeV. The 
chiral extrapolated results are G^^(O) = —0.017(25), 
Aa = 0.58(16), (r^)^./ = -7.4(71) x lO^^fm^ and g% = 
0.027(16) (or (r^)^; = -2.4(15) x lO'^fm^). 

Before quoting the final results, we consider the sys- 
tematic uncertainties yet to be addressed. First, we 
analyze the ambiguity of dependence in form fac- 
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FIG. 3: The chiral extrapolated results for G%}{Q'^) (upper) 
and Ge{Q^) (lower) plotted with solid lines. Shaded regions 
represent the error-band with statistical and systematic error 
added in quadrature. Shown together are the lattice data (and 
Q^-extrapolated Gm{0)) for Kud = 0.13760 (circles), 0.13800 
(triangles), 0.13825 (squares) with offset for visibility. 



tors, which is rather unknown for the strangeness. We 
also study the monopole form, G|^(Q^) = G'|^(0)/(1 -f 
QVA2), G%{Q^) = g% ■ 0V(1 + <3VA'), and find that 
the precision of the lattice data cannot disentangle the 
difference of the dipole/monopole behaviors. The results 
from the monopole fit are G|^(0) = -0.024(39), Aa = 
0.34(17), {rl)Ai = -12(18)xl0-3fm^ and g|j = 0.035(17) 
(or (r^)B = —3.1(15) x 10^"^fm^), which are consistent 
with those from the dipole fit and have relatively larger 
statistical errors. 

Second, we study the uncertainties in chiral extrapola- 
tion by testing two alternative extrapolations. In the first 
one, we take into account the nucleon mass dependence 
on the quark mass, using the lattice nucleon mass. From 
the physical viewpoint, this corresponds to measuring the 
magnetic moment not in units of lattice magneton but 
physical magneton [H]. We obtain G^,^(0) = -0.017(24), 
Aa = 0.77(21), {r'^^)M = -4.6(43) x lO-^fm^ from the 
dipole fit, and G|,f(0) = -0.023(36), Aa = 0.45(22), 
{tDm = ^7(11) X lO^^fm^ from the monopole fit, while 



G|; cannot be fitted because fit dof are not sufficient. 
Note that results are consistent with previous analy- 
ses. In the second alternative, we use the linear fit in 
terms of m^, observing the results have weak quark 



{ri)M = -3(17) X 10-^fm^ g% = 0.023(12) (or (rf)^^ = 
-2.0(10) X IQ-Hm^) from the dipole fit, and G^(0) = 
-0.023(33), {t'Dm = -8(43) x lO-^fm^, g% = 0.029(12) 
(or {rfjE = —2.5(11) x lO^'^fm^) from the monopole fit, 
while pole masses are found to suffer from too large statis- 
tical errors to extract useful information. We find again 
the results are consistent with previous analyses. While 
a further clarification with physically light quark mass 
simulation and a check on convergence of HBxPT [26] is 
desirable, we use the dependence of results on different 
extrapolations as systematic uncertainties in the chiral 
extrapolation. 

Third, we examine the contamination from excited 
states. Because our spectroscopy study indicates that the 
mass of Roper resonance is massive compared to the 5*11 
state on the current lattice (27| . the dominant contami- 
nations are (transition) form factors associated with 5'ii. 
On this point, we find that such contaminations can be 
eliminated theoretically, by making the following substi- 
tutions: r± ^ f± = (1 ± mN'/E%, ■ 74) /2 in Eq. ©, 



r- 



Here, totv* denotes the 5*11 mass, 
f^" = {niN* T *;^-)/(2tojv*) with p'_ 



in Eqs. (HI) and ©. 



{E%.,p'). Note 
also that the (r.h.s.) of Eqs. ([4]) and (O have modifi- 
cations in the kinematical factor, which we do not show 
explicitly. It is found that the results from these equa- 
tions are basically the same as before, so we conclude that 
the contamination regarding the Sn state is negligible. 

As remaining sources of systematic error, one might 
worry that the finite volume artifact could be substan- 
tial considering that the physical spacial size of the lattice 
is about (2fm)^. However, we recall that Sachs radii are 
found to be quite small, \{r'^)E,M\ <C O.lfm^, which indi- 
cates a small finite volume artifact. For the discretization 
error, we first examine the error associated with the fi- 
nite momentum q, using the dispersion relation of the 
nucleon. We find that the nucleon energy at each q on 
the lattice is consistent with the dispersion relation, and 
conclude that finite (qa) discretization error is negligible. 
As another discretization error, we note that rriN {itik) 
is found to have 6 (8) % error for the current configura- 
tions [1^, [2^ . Considering the dependence of G|; on 
these masses, we estimate that the discretization errors 
in our results amount to ^ 10%, and are much smaller 
than the statistical errors. Of course, more quantitative 
investigations on these issues are necessary with larger 
and finer lattices, and such work is in progress. 

Here, we present our final results. For the magnetic 
moment, Gi^(O) = -0.017(25)(07), where the first er- 
ror is statistical and the second is systematic from un- 
certainties of the extrapolation and chiral extrapola- 
tion. We also obtain Aa — 0.58(16) (19) for dipole mass 
or Aa — 0.34(17) (11) for monopole mass, and g|; = 
0.027(16)(08). These lead to GljiQ"^) = -0.015(23), 
G%{Q^) = 0.0022(19) at = 0.1GeV^ where error is 
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obtained by quadrature from statistical and systematic 
errors. Note that these are consistent with the world av- 
eraged data 0, H, 01 , with an order of magnitude smaller 
error. In Fig. [3 we plot G%j{Q^), G%{Q^), where the 
shaded regions correspond to the square-summed error. 
The lattice data for each Kud are also plotted. 

IV. SUMMARY 

We have studied the strangeness electromagnetic form 
factors of the nucleon using Nf = 2 + 1 clover fermion 
configurations. In the evaluation of the disconnected in- 
sertion (DI), the Z(4) stochastic method along with unbi- 
ased subtractions from the hopping parameter expansion 
has been used. We have developed several techniques to 
achieve a good S /N in the DI calculation, and found that 
increasing the number of nucleon sources for each config- 
uration is particularly efficient. For all quark mass simu- 
lations, G%[{Q^) has been found to be negative with 2-3 a 



signals for low regions. Upon and chiral extrapola- 
tions, we have obtained GljiQ"^ = 0) = -0.017(25) (07), 
where the first error is statistical, and the second reflects 
the uncertainties in and chiral extrapolations. We 
have also obtained G'm(Q^) = -0.015(23), G^Q^) = 
0.0022(19) at O.lGeV^, where error is obtained by 

quadrature from statistical and systematic errors. These 
results are consistent with experimental values, and our 
errors are an order of magnitude smaller. 
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